Electrical neutrality and pion modes in the two flavor PNJL model 



H. Abuki,i'E|M. Ciminale,i'2.[t|R.Gatto,3.[lN. D. Ippolito,!- 2.§ G. NarduUi,!- and M. Ruggierii. 2.0 

'I.N.F.N., Sezione di Bam, 1-70126 Ban, Italia 
^Universita di Bari, 1-70126 Bari, Italia 
^ Departement de Physique Theorique, Universite de Geneve, CH-1211 Geneve J,., Suisse 

We study the phase diagram and the pion modes in the electrically neutral two flavor PNJL 
model. One of the main result of this paper is that in the model with massive quarks, when 
electrical neutrality is required, pious do not condense in the ground state of the model: the isospin 
chemical potential /i/ = —fie/ 2 is always smaller than the value required for pion condensation to 
occur. Moreover we investigate on the pious and a mass spectra. We find that the qualitative 
behavior of the masses resembles that obtained in the NJL model. We close this paper by studying 
the intriguing possibility that bound states with the quantum numbers of the pious can be formed 
above the chiral phase transition. 



I. INTRODUCTION 

The use of efFective models to understand the phases of Quantum Chromodynamics (QCD) is nowadays a popular 
tool. Among them, the Nambu-Jona Lasinio model (NJL) is widely used since it allows for a simple discussion of 
chiral symmetry breaking [2]. The main defect of the NJL model is that it completely ignores glue dynamics: 
gluons are not introduced in the lagrangian, and the QCD quark-quark interaction mediated by gluons is replaced by 
an efFective four fermion interaction. 

An improvement of the NJL model is provided by the introduction of a background temporal static gluon field, 
coupled to quarks via the QCD covariant derivative. The gluon field is related to the expectation value of Polyakov l oop 
$ d, and the new model is called Polyakov-NJL (PNJL in the following) 0, i, i, 0, i, i, 113, [HI, 113, [11 Q 111, [M lizl 
It is well known that the Polyakov loop serves as an order parameter for confinement deconfinement transition in 
the pure glue theory In more detail, the confined phase is characterized by ($) = and a Z3 symmetry; on the 
other hand in the deconfined phase ($) ^ and the symmetry Z3 is broken. In presence of dynamical quarks, Z3 
symmetry is explicitly broken and $ can not be used as an order parameter. Nevertheless it is commonly used as an 
indicator of the deconfinement transition. 

In the PNJL model an effective potential for $ is added by hand to the quark lagrangian, and $ is coupled to the 
quarks via the QCD covariant derivative. The value of ($) in the ground state, as well as other quantities of interest 
like the constituent quark masses, are obtained by minimization of the thermodynamic potential, or equivalently by 
solving coupled gap equations as in the NJL case. Despite the apparent complication due to the increase of the 
degrees of freedom, the PNJL model offers a better description of QCD then NJL does since it allows to derive, in 
the framework of field theory, several results obtained in lattice simulations [3, 0, H, [3, [HI [13, [3 • 

The knowledge of the phase diagram of the PNJL model as well as of its bulk properties is of a certain interest. In 
this paper we study the phase diagram and the pion modes in the electrically neutral two flavor PNJL model. The 
motivation of our study is straightforward: we are interested to investigate the equilibrium ground state of the model, 
and in the ground state a net electric charge can not be present, whichever the value of the temperature and of the 
quark chemical potential is. We consider the possibility of chiral as well as pion condensation in the ground state, 
at zero and non zero quark chemical potential. A related study without electrical neutrality and at zero chemical 
potential has been done in Ref. 14]; see also Ref. 18] for a similar study in the framework of the electrically neutral 
NJL model in the chiral limit. 

One of the results of our work is that when electrical neutrality is required, pion do not condense in the ground 
state of the model: the isospin chemical potential /i/ = — /Xe/2, with /ig the electron chemical potential, is always 
smaller than the value required for pion condensation to occur. Moreover we investigate on the pions and a mass 
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spectra. We find that the quahtative behavior of the masses resembles that obtained in the NJL model. We close this 
paper by studying the possibility that a bound state with the quantum numbers of the pions can be formed above 
the chiral phase transition. 

The plan of the paper is as follows: in Section [IT] wc briefly describe the PNJL model; in Section Hill we show the 
phase diagram of the model, and some results for quantities of interest (in medium quark mass, pion condensate, mean 
value of the Polyakov loop, electron chemical potential). In Section lTVl we compute the masses of the pseudo-Goldstone 
modes of the model and of the a mode. Finally in Section |V] we draw our conclusions. 

II. FORMALISM 

The Lagrangian density of the two flavor PNJL model is given by 

C = e{i-i^d^)e + ^ {i-i^D^ + /170 - m) V + G [(i^V)' + {4'^l^T^pf'\ - U[^, $, T] . (1) 

In the above equation e denotes the electron field; ip is the quark spinor with Dirac, color and flavor indices (implicitly 
summed), m corresponds to the bare quark mass matrix; we assume from the very beginning rriu = rrid- The 
covariant derivative is defined as usual as = — iAfj_. The gluon background field A^^ = 6q^Aq is supposed 
to be homogeneous and static, with ^0 = g^oTa and Ta, a = 1, ... ,8 being the SU{3) color generators with the 
normalization condition Tr[rarb] = Sab- t is a vector of Pauli matrices in flavor space. Finally fi is the chemical mean 
quark chemical potential, related to the conserved baryon number. 

In Eq. ([1]) $, $ correspond to the normalized traced Polyakov loop L and its hermitian conjugate respectively, 
$ = TrL/TVc, i = TiLyNc, with 

L = 'Pcxp^ J Aidrj = exp (i/SAi) , Ai = iAo , (2) 

and /3 = 1/T. The term U[^,^,T] is the effective potential for the traced Polyakov loop; in absence of dynamical 
quarks it is built in order to reproduce the pure glue lattice data of QCD, namely thermodynamical quantities 
(pressure, entropy and energy density) and the deconfinement temperature of heavy (non-dynamical) quarks, T = 270 
MeV. Several forms of this potential have been suggested in the hterature, see for example ^, iiS]- In this paper 
we adopt the following logarithmic form Q, 

T] ^ _b2(p_^^ ^ ^^y^ _ ^ ^^1,3 _^ ^3) „ 3(|,$)2] ^ (3) 

with 

62(r)=ao + ai(^^)+a2(^^) , 5(T) = 63 (y) ' (4) 

Numerical values of the coefficients are as follows [3] : 

ao = 3.51 , ai = -2.47 , as = 15.2 , 63 = -1.75 . (5) 

If dynamical quarks were not present then one should chose Tq = 270 MeV in order to reproduce the deconfinement 
transition at T = 270 MeV. In p resence of quarks Tq gets a dependence on the number of active flavors, as shown in 
Refs. d, [l3l- Following Ref. |l0l | we chose Tq = 208 MeV in Eq. (|4]), which is appropriate to deal with two degenerate 
flavors. 

As explained in the Introduction we are interested to study the ground state of the model specified by the Lagrangian 
in Eq. (HI), at each value of the temperature T and the chemical potential /x, corresponding to a vanishing total electric 
charge. To this end we introduce electrons in Eq. ([1]) since a net number of electrons could be needed, beside the 
pion condensate (if any, see below), to ensure the electrical neutrality of the ground state. In order to build the 
neutral ground state we work as usual in the gran canonical ensemble formalism, adding to Eq. ([1]) the term iiqNq, 
fiQ being the chemical potential for the total charge Nq, and requiring stationarity of the thermodynamic potential 
with respect to variations of fig, which is equivalent to the requirement < Nq >= in the ground state. To be more 
specific, since the total charge operator is given in terms of the quark and electron fields by 

Nq = '^Nu - ^Nd - iVe = ^u^u ~ ^dU - ete , (6) 
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it is easy to recognize that the lagrangian £ in the gran canonical ensemble C = C + iiqNq can be written as 

£ = e(i7^9^ + /ie7o)e + ^ {i-ft^D^ + fi-fo - m) + G (i/^V')^ + (V^^TsrV')^ - ^['^>, ^, T] , (7) 
where /ie = —^iQ and the quark chemical potential matrix fl is defined in flavor-color space as 

.. ,V, "l^le, (8) 



^+3^e 

where Ic denotes identity matrix in color space. 

In this paper we work in the mean field approximation. In order to study chiral symmetry breaking and to allow 
for pion condensation we assume that in the ground state the expectation values for the following operators may 
develop [Hill, 

a = (Tp'ip) , TT = (i/Ji75riV') . (9) 

In the above equation a summation over flavor and color is understood. We have assumed that the pion condensate 
aligns along the ti direction in flavor space. This choice is not restrictive. As a matter of fact we should allow for 
independent condensation both in tt"*" and in tt^ channels 

7r+ EE (V;i75T+V^) = -^e'^ , TT- = {i;i-f5T.i}) = ^e"** , (10) 

with T± — (ti ± T2)/^/2; but the thermodynamical potential potential is not dependent on the phase 6, therefore we 
can assume 9 = which leaves us with 7r+ = tt^ = 7r/V2 and introduce only one condensate, specified in Eq. ([9|). In 
the mean field approximation the PNJL lagrangian reads 

£ = e(^7^a^ + fielo)e + ^ (^7^^'' + ^70 - M - Ni-fsTi) - G [a^ + tt^] - U[^, !>, T] , (11) 

where the in-medium quark mass M and the pion field N are related to a and tt by means of the relations 

M = m- 2Ga , (12) 
N = -~2Gtt . (13) 

The thermodynamical potential f2 can be obtained by integration over the fermion fields in the partition function of 
the model, see for example Ref. (Toj . 

where the sum is over fermion Matsubara frequencies ujn = TTT(2n + 1), and the trace is over Dirac, fiavor and color 
indices. A is an ultraviolet cutoff to ensure convergence of the momentum integral. The inverse quark propagator is 
defined as 

5-^(-«,p)=f + + ^ ^1 AfVlc (15) 

The ground state of the model is defined by the values of cr, tt, $, $ that minimize f2 and that have a vanishing total 
charge; the latter condition is equivalent to the requirement 

f = 0. (16, 

The parameters m, G and A are given by 0] 

m = 5.5 MeV , G = 5.04 GeV"^ , A = 650.9 MeV , (17) 

which fix, at zero temperature and zero chemical potential, the pion mass ttZtt = 139.3 MeV, the pion decay constant 
/tt = 92.3 MeV and the chiral condensate (uu) = -(251 MeV)^. 
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FIG. 1: Phase diagram of the electricaUy neutral two flavor PNJL modeL Dashed hne corresponds to the chiral crossover; 
sohd hne describes the first order chiral transition. At ^ = the critical point is located at T = 206 MeV; at T = the chiral 
transition is found at = 348 MeV. The black dot denotes the critical end point (CEP), located at {he,Te) — (342, 79) MeV. 




FIG. 2: Comparison of the first order transitions of the PNJL model with neutrality (called Model I in the text), of the PNJL 
model with /j,e = (model II) and of the NJL model with electrical neutrality (model III). 



III. PHASE DIAGRAM OF THE MODEL 



In Fig. [T] we plot the phase diagram of the electrically neutral two flavor PNJL model. Dashed line corresponds to 
the chiral crossover; solid line describes the first order chiral transition. At each value of fi the crossover is identified 
with the inflection point of the chiral condensate. Analogously the first order transition is defined by the discontinuity 
of cr. In the region below the lines the chiral symmetry is broken, above the lines it is restored. At /x = the chiral 
symmetry is restored aXT — 206 MeV. For comparison, the inflection point of the Polyakov loop (which is commonly 
associated to the deconfinement crossover) at ^ = is located at T = 180 MeV. Moreover, at T = the chiral phase 
transition is of first order and is found at /j, = 348 MeV. The black dot denotes the critical end point (CEP), located 
at {fiE,TE) = (342,79) MeV. As explained later, we find a vanishing pion condensate once electrical neutrality is 
required. 

It is instructive to compare the PNJL results discussed above with the electrically neutral NJL model, as well as the 
PNJL model with /ig = 0, in order to emphasize the role of electrical neutrality and of the Polyakov loop. To this end 
in Fig. [2] we compare the first order transitions of the PNJL model with neutrality (model I), of the PNJL model with 
/Xe = (model II) and of the NJL model with electrical neutrality (model III). Comparison of models I and II shows 



5 



(1=0 MeV (1=342 MeV 





1.0 
0.8 






1.0 
0.8 












ivr\ 






M 










0.6 




z" 


0.6 














0.4 






0.4 














0.2 
0.0 






0.2 
0.0 




N 















50 100 150 200 250 50 100 150 200 250 

T(MeV) T(MeV) 



FIG. 3: Quark constituent mass A4 (in units of Mo), pion condensate A'^ (in units of Mo) and $ as a function of temperature 
(in MeV) for different values of fi. Mo = 325 MeV is the constituent quark mass evaluated at ^ = 0, T = and /Xe = 0. 

that adding the electron chemical potential and requiring electrical neutrality enlarges of some MeV the /x window of 
the chirally broken phase. Moreover, comparison of models I and III shows that adding the self-consistently evaluated 
Polyakov loop to the model renders the broken phase more robust, as the critical temperatures are increased. 

In Fig. [3] we plot the quark constituent mass M (in units of Mq), pion condensate N (in units of Mq) and <& as 
a function of temperature (in MeV) for two values of /i. We have defined Mq — 325 MeV as the constituent quark 
mass evaluated at /Lt = 0, T = and He = 0. The electron chemical potential evaluated at some values of ^ is shown 
in Fig. m We observe that when electrical neutrality is required, pions do not condense in the ground state of the 
model. As a matter of fact TV = in Fig. [Sj we have verified this result for several values of /i. Stated in other words 
the isospin chemical potential fj,i = in the neutral phase is always smaller than the value required for pion 

condensation to occur. 

We have verified the robustness of our calculations in several ways. First, we have treated as a free parameter, 
reproducing the results of Ref. [14i] finding a pion condensate for high values of /ig. Second, we have compared 
our results with those Ref. [l^. The authors of Ref. [l8| study the neutral ground state in the NJL model at low 
temperature. They use a set of parameters in which the values of G and A are similar to ours, but they consider only 
the chiral limit m„ = nid = 0. In this limit they find a narrow window in ^ at zero temperature in which the pion 
condensate can exist. We have explicitly checked that putting by hand in-u = rud = in our calculations we reproduce 
the aforementioned window. This is a further check since at low temperature the PNJL reduces to the NJL model. 
Moreover we have verified that the pion condensate is washed out if the current quark mass is increased from zero 
to the physical value r7i„ = nid — 5.5 MeV. This is explicitly shown in Fig [5] where we plot the pion condensate N 
and the constituent quark mass M in units of Mq, the constituent quark mass at T = ^ = 0, as a function of the 
bare quark mass tuq at the representative value of /i = 330 MeV (we have checked that the same result hold for other 
values of /i) and T = 0. Therefore the current quark mass works as an external field that drives the alignment of the 
vacuum along the chiral condensate direction. We stress that Fig. 5 has to be taken only as a trick that allows a clear 
comparison between our results and those of Ref. [ll] . A true plot of M and as a function of the current quark 
mass requires the self consistent evaluation of G and A, which is beyond the scope of the present paper since we study 
only the phase diagram and the pion modes at the physical point. Nevertheless the problem of the modification of 
the ground state by variations of the current quark mass is interesting and deserves further study, therefore we leave 
this to a future project. 

IV. PSEUDO-GOLDSTONE MODES 

In this Section we compute the masses of the pseudo-Goldstone modes as well as of the a mode in the two flavor 
and electrically neutral PNJL model. As discussed in the previous Section for each value of fi and T, after requiring 
electrical neutrality, the difference of chemical potentials between up and down quarks is lower than the critical value 
required for the realization of the pion condensation. In this Section we therefore put by hand A^ = in the quark 
propagator, Eq. (fT51) . 

The pion mass is defined as the solution of the well-known Bethe-Salpeter equation for the bound states. Its 
derivation in the context of the NJL model is in Ref. 0. The same equation is valid in the PNJL model, the 
difference between PNJL and NJL being only the different quark propagator to be used in the calculation of the pion 
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FIG. 4: Electron chemical potential as a function of the temperature, for different values of the quark chemical potential. 
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FIG. 5: Pion condensate A'^ and the constituent quark mass M in units of Mo, the constituent quark mass ai T — jj, — 0, as & 
function of the bare quark mass mo at the representative value of = 330 MeV and T = 0. 



polarization tensor. The equation reads 



1 ~ 2GIlA{m\) ^ , A = 1,2, 3 



(18) 



In the above equation 11^ is the polarization tensor of the pion A, specified later, where A = 1,2,3 correspond 
respectively to tt^, tt" and tt''; itia denotes the pion mass, and G is the coupling constant introduced in Eq. ([7|. 
The solutions of Eq. ((TH]) correspond to the poles of the pion propagator, the latter evaluated in the random phase 
approximation 

We now specify the pion polarization tensor. We introduce the polarization matrix 



n,,(m2) 



Tr 



(19) 



here i,i = 1,2,3, — (m, 0), and the trace is over color, flavor and Dirac indices; the matrices fi are operators in 
the flavor space, defined by 



Tl,2 



1 



(n ± iT2) , Ta = Ta , 



(20) 



where f^, i = 1,2,3 are the Pauli matrices; S is the fermion propagator, whose inverse is defined in Eq. ([T51) with 
N = 0, and masses and chemical potentials evaluated self consistently in the electrically neutral phase. In terms of 
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FIG. 6: Masses of the charged pions (left panel) and of the neutral pion (right panel) in the electrically neutral phase, as a 
function of the temperature, for different values of the quark chemical potential. 




FIG. 7: Mass of the a mode in the electrically neutral phase, as a function of the temperature, for different values of the quark 
chemical potential. 



the matrix Ily the pion polarization tensors appearing in Eq. (jl9p are defined as 



Hi = ni2 



n. 



33 



Analogously the equation for the a mass reads 

l-2Gn,(m2) = 



with the polarization tensor defined as 



n,(m2) 



Tr 



iS p 



(21) 



(22) 



(23) 



Once the traces in Eqs. and are evaluated, we are left with the loop 4-momentum integral. We replace 
the integral over energy by a summation over fermion Matsubara frequencies; moreover, the integral over directions 
of momentum is trivial. Therefore we are left with the integration over the modulus of p only, which is evaluated 
numerically. 

In Fig. [6] we plot the masses of the charged pions (left panel) and of the neutral pion (right panel) in the electrically 
neutral phase, as a function of the temperature, for different values of the quark chemical potential. The behavior 
of the pseudoscalar modes as the temperature is increased is qualitatively the same observed in the NJL model 0] ■ 
At each value of fi, the charged and neutral pion masses are of the order of the zero temperature value « 140 MeV 
below the chiral crossover (or the chiral first order transition); as the chiral transition (either first order or crossover) 
occurs, the pion masses rapidly increase signaling the disappearing of the pseudoscalar modes from the low energy 
spectrum of the model. 

In Fig. [7] we plot the mass of the a mode in the electrically neutral phase, as a function of the temperature, for 
different values of the quark chemical potential. As in the case of the pions, the behavior of m„ is analogous to that 
observed in the NJL model 
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FIG. 8: The ratio m-^/2M, with Al the constituent quark mass, as a function of temperature, for two representative values of 
the quark chemical potential. The vertical bold dashed line denotes the chiral transition temperature. 
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TABLE I: Comparison of the pion melting temperatures Tm and the chiral transition temperatures T,,, with electrical neutrality 
and with /le = 0. Quark chemical potentials and temperatures are measured in MeV. 

It is interesting to evaluate the ratio m7r/2M, with M the constituent quark mass, as a function of temperature. 
If it is larger than one, than a bound state with the quantum numbers of the pion is less stable than a state made of 
a free quark and a free antiquark, and thus the pion melts to its constituent quarks. In Fig. [8] we show the results 
of the computation of mT,/2M at different values of /i in the electrically neutral phase, for values of temperatures 
close either to the chiral crossover or to the first order chiral transition. We notice that for any value of the quark 
chemical potential there exists a window in temperature, above the chiral restoration temperature, in which the pions 
are lighter than a quark and an antiquark. The window shrinks as /i is increased. These results show that even above 
the critical temperature a bound state with the pions quantum numbers can be formed. 

We observe that the window in which the bound state is stable above the chiral phase transition (as well as the 
chiral crossover) is not a, peculiarity of the PNJL model. As a matter of fact the same window was noticed in the NJL 
model at /i = in Ref. [23l- We have verified that it exists in the neutral NJL model as well. For example, at /.t = 
we find that the bound states exist in the window (195, 212) MeV; this result has to be compared with the PNJL one, 
namely (205, 245) MeV, showing that the PNJL window is almost 2.4 times larger than the NJL one. Analogously at 
fi = 342 MeV the NJL window is (34, 44) MeV, to be compared with the PNJL result (73, 85) MeV. 

Similar results have been obtained in the framework of the PNJL model in Ref. [21|. The authors of Ref. [2l| 
consider the two flavor PNJL model without the constraint of electrical neutrality. Moreover they use a polynomial 
form of the Polyakov loop effective potential U in Eq. ([7]) instead of the logarithmic one, see Eq. Therefore a 
quantitative comparison with their results is not feasible. Instead we have computed the ratio Tm/T^ in the case 
/ie = in the PNJL model with the logarithmic form of the effective potential, for several values of in order to 
quantify the role of the electrical neutrality on the ratio Tm /T-^. Our results are summarized in Table HI 

It is interesting to superimpose the bound state existence region to the phase diagram in Fig. [T] This is done 
in Fig. [9] where we show the region in the /i — T plane where the bound state can be formed. In the white region 
below the gray domain the chiral symmetry is broken and pions live as pseudo-Goldstone modes. In the gray region 
chiral symmetry is restored but mTr/2Afg < 1, thus bound states can be formed. Finally, in the white region above 
the gray domain chiral symmetry is restored and free quark states are more stable than bound states (this region is 
characterized by 771,^/2 Af^ > !)• 

The existence of bound states above the chiral transition temperature is interesting because it implies that quark 
matter above the chiral transition is strongly coupled. It has been suggested [l^, [Ij] that the formation of bound 
states just above the chiral transition can explain the experimental results obtained by non-central heavy ion collisions 
at RHIC facility [l^, H^, 11^ [l^] ; in particular it provides a simple explanation to the low viscosity-to-entropy ratio 
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observed in elliptic flow simulations. The next step would be the computation of the transport coefficients of quark 
matter just above the chiral transition temperature. We leave this to a future project. 
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FIG. 9: Region of existence of bound states in the i-i — T plane. x^B denotes the region where chiral symmetry is broken; 
xSR denotes the region where chiral symmetry is restored. M is the constituent quark mass. The dashed and the solid lines 
correspond respectively to the chiral crossover and to the chiral first order transition, the dot denotes the CEP. In the gray 
region chiral symmetry is restored but a bound state with the quantum numbers of the pious can still be formed. See the text 
for more details. 



V. CONCLUSIONS 



In this paper we have studied the two flavor PNJL model in presence of a charge chemical potential, requiring 
electrical neutrality in the mean field approximation. We have assumed that in the ground state both the chiral and 
the pion condensates develop, and we have evaluated their values, as well as the value of the Polyakov loop, for each 
value of the temperature and of the quark chemical potential by minimizing the thermodynamic potential under the 
condition that the total electric charge of the system vanishes. 

We have drawn the phase diagram of the model, see Fig. [TJ One of our main results is that we do not observe pion 
condensation in the mean field approximation. We relate this to the non vanishing current quark mass, as shown in 
Fig. 5. This part of our study is completed by the results for the constituent quark masses, the Polyakov loop and 
the charge chemical potential, see Figs. [3] and S) 

We have computed the masses of the pseudo-Goldstone modes and of the cr-mode. The behavior of the masses as a 
function of the temperature for different values of the quark chemical potential are shown in Figs.[6]and[7l Furthermore 
we have investigated on the possibility of existence of bound states with the quantum numbers of the pions above 
the chiral critical temperature. To this end we have compared the computed pion mass with twice the constituent 
quark mass. The result is shown in Fig. [8] We have found that for any value of the quark chemical potential there 
exists a window in temperature, above the chiral restoration temperature, in which the pions are lighter than a state 
formed by a free quark and a free antiquark. The aforementioned window shrinks as fj, is increased. Our results show 
that even above the critical temperature a bound state with the pions quantum numbers can be formed. This is 
summarized in Fig. [9l 

In this work we have not studied independently the role of the chemical potential on the suppression of the pion 
condensate. This can be done by leaving the isospin chemical potential as a fixed and free parameter, and varying 
only the mean quark chemical potential. A further natural extension of our work is the study of the 2 + 1 flavor 
PNJL model with meson condensation, with and without neutrality conditions implemented. Moreover, it would be 
interesting to include a color superconductive state in the neutral model [1^ (studies of the superconductive state 
without neutrality conditions can be found in Refs. 

U,M)- Even more one should compute the transport coefficients 
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of quark matter just above the chiral transition temperature. We leave these projects to future works. 
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